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ABSTRACT 

Flows of synchrotron-emitting material can be found in several astrophysical contexts, including extragalac- 
tic jets and pulsar-wind nebulae (PWNe). For X-ray synchrotron emission, flow times are often longer than 
1 electron radiative lifetimes, so the effective source size at a given X-ray energy is the distance electrons radiat- 

| ing at that energy can convect before they burn off. Since synchrotron losses vary strongly with electron energy, 

the source size drops with increasing X-ray energy, resulting in a steepening of the synchrotron spectrum. For 
homogeneous sources, this burnoff produces the well-known result of a steepening by 0.5 in the source's in- 
tegrated spectral index. However, for inhomogeneous sources, different amounts of steepening are possible. 
I exhibit a simple phenomenological picture of an outflow, with transverse flow-tube radius, magnetic-field 
^vq ■ strength, matter density, and flow velocity all varying as different powers of distance from the injection point. 

For such a picture, I calculate the value of the spectral index above the break as a function of the power-law 
indices, and show the possible range of steepenings. I show that these simple calculations are confirmed by 
full integrations of source luminosity, which also include the spectral "bump" below the break from the accu- 
mulation of electrons formerly at higher energies. In many cases, extragalactic jets show X-ray synchrotron 
emission steeper by more than 0.5 than the radio emission; the same phenomenon is exhibited by many pulsar- 
wind nebulae. It is possible that source inhomogeneities are responsible in at least some cases, so that the 
amount of spectral steepening becomes a diagnostic for source dynamical or geometrical properties. 
Subject headings: galaxies: jets — radiation mechanisms: non-thermal — supernova remnants — supernova 
remnants: individual (B0540-693) — X-rays: ISM 

a 

1. INTRODUCTION 

1.1. Spectral breaks in jets and pulsar-wind nebulae 

Broad-band spectra from synchrotron radiation characterize a wide variety of astrophysical sources, including active galactic 
I/"") , nuclei (AGN), shell supernova remnants (SNRs), and pulsar-wind nebulae (PWNe). When observed over a sufficiently broad 
C — ■ ■ frequency range, such spectra almost invariably show steepening to higher energies. Such steepening can be attributed either to 
intrinsic spectral structure or to the effects of radiative losses. However, the power-law shape of radio spectra {S u oc v~ a ), with 
, a ~ 0.5-0.8 for many sources, suggests an origin of the requisite relativistic electrons in diffusive shock acceleration (DSA), 
which produces a power-law (or near-power-law, for efficient nonlinear DSA) spectrum of particles N(E) oc E~ s with s = 2a+ 1 
° depending on the shock compression ratio. Shocks are clearly present in these astrophysical objects: outer blast waves in SNRs 
C"^ . (and perhaps reverse shocks into ejecta as well, for younger objects); relativistic -wind termination shocks in PWNe; and shocks 
^ ■ in jets and hotspots in active galaxies. The absence of an obvious mechanism to generate a broken power-law distribution, and 
the necessity of the operation of radiative losses at some level, has led to the common acceptance of synchrotron losses as the 
mechanism to bring about spectral steepening. 

The standard calculation of the effect of synchrotron losses (reviewed in the next section) for a homogeneous source predicts 
a steepening of the initial electron spectrum to a second power-law one power steeper (*2 = s+1) than the injection spectrum, 
implying a radiation spectrum one-half power steeper (at = a + 0.5). This is in fact rarely observed. In shell supernovae, the 
maximum energies to which electrons can be accelerated are limited by losses or by finite shock age (or size), but produce 
an exponential cutoff in N(E), further broadened by inhomogeneities, as observed in a few cases in which synchrotron X-ray 
emission can be identified (see Reynolds 2008 for a review), so a sharp spectral break to a steeper power-law is neither expected 
nor observed. In PWNe, radio spectra are almost all flatter than a = 0.5, a phenomenon not well explained at present, but the 
steepenings A = ahigh _ aiow are almost always greater than 0.5 (ranging from 0.7 to 1, in seven of the eight cases tabulated in 
Chevalier 2005, using updated values for B0540-693 from Williams et al. 2008). Knots in extragalactic jets, when observable 
in X-rays, show similar too-large steepenings (e.g., M87: Perlman & Wilson 2005; Cygnus A: Stawarz et al. 2007). One 
standard interpretati on of knot and hot-spot spectr a invokes shock acceleration and subsequent spectral steepening in a uniform 
post-shock region dHeavens & Mei senheimer 1987]), but it cannot explain these larger values of A. Coleman & Bicknell (1988) 
present numerical hydrodynamic simulations and find larger values of A, which they apply to observations, but without analytic 
results allowing the wider application of the results. While there is recognition of the possibility of values of A different from 
0.5 (e.g., Kennel & Coroniti 1984a, who find A = (4 + a)/9 for the Crab Nebula, or Petre et al. 2007, in a qualitative discussion 
of the broadband spectrum of the PWN B0540-693), there has as yet been no simple characterization of conditions under which 
values of A y 0.5 can naturally arise. That characterization is the goal of this paper. 
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1.2. Synchrotron losses 
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The first widely known calculation of the behavior of a distribution of electrons subject to synchrotron losses is that of Karda- 
shev (1962). While most of these results are well known, it is important to recall the particular conditions under which each is 
applicable, so I shall beg the reader's indulgence for a brief review, which can also serve to fix notation. Kardashev solved the 
kinetic equation for a distribution N(E,t) of electrons subject to gains by first and second-order Fermi acceleration, and losses 
due to radiation or adiabatic expansion, with and without the assumption of new-particle injection and stationarity. He writes the 
energy-loss rate from a single electron as 

E = -bB 2 L E 2 , (1) 

where B± = Bsind, b = (2/3)(e 4 /m 4 c 7 ) = 2.37 x 10~ 3 cgs (e.g., Pacho lczykll 1 970h . and 9 is the electron pitch angle between 
its velocity vector and the magnetic field. Kardashev pointed out that in a uniform magnetic field in the absence of scattering, 
electrons preserve their value of 9, since they radiate a beam pattern that is symmetric with respect to their velocity vector. An 
electron injected into B at t = with energy Eq has an energy E after time t given by the well-known result 

E(t)= . . (2) 

An initially infinitely energetic electron is reduced after time t to energy E max (t,9) = 1 /bB 2 ± t. A power-law energy distribution of 
electrons N(Eq) = KEq s , with a single value of pitch angle, will be cut off at E m - dx (t,9). The electron distribution N(E) will evolve 
according to N(E)dE = N(Eo)dEo, so that 

N(E) = K[E (E)r^, (3) 
dE 

with E(Eo) given by Equation|2] If s < 2, the electrons initially above E mlix {t) are sufficiently numerous to pile up in a "bump" just 
below £m ax (fX while if s > 2, the "bump" disappears. (Numerical solutions to Equation [3] are shown below for inhomogeneous 
models, illustrating the "bump" effect.) 

An initially isotropic distribution of electrons suffers unequal radiation losses, with electrons with large pitch angles being 
more rapidly depleted. For an initial isotropic power-law distribution, after a time t one finds no electrons with pitch angles 
greater than 6* max given by sin 2 # max = 1 / (bEB 2 t). Since for synchrotron radiation, an individual electron's radiation pattern has an 
angular width of order 1 where 7 is the individual Lorentz factor, and 7 > 10 3 for radio emission and higher frequencies, we 
can approximate electrons as radiating exactly in their directions of motion. Then 9 is also the angle between the line of sight and 
the local magnetic field. An initially isotropic distribution of electrons injected at t = into a source with a uniform magnetic field 
making an angle \ with the line of sight, observed through its synchrotron emission, would then disappear abruptly at time t(~x). 
More realistically, one might expect that the source has a tangled magnetic field, so that all values of \ are achieved in some part 
of the source. One should then (for an unresolved source) perform an angular integration over the electron distribution. The result 
is an electron distribution that steepens by one power of E, i.e., N(E) oc E~"~ l , above the characteristic energy Et,= l /(aB 2 t). (The 
synchrotron emission from this distribution steepens above v{Et,) by more than the value of 0.5 in spectral index a (S v oc v~ a ) 
because, in this time-dependent case without continuous injection, a correlation exists between E and 9 such that more efficiently 
radiating electrons are depleted most rapidly.) 

This situation is still relatively unrealistic, as it ignores any processes by which electrons could change their pitch angles. If 
electrons scatter in pitch angle on timescales much shorter than the synchrotron-loss timescale, one should simply average the 
energy-loss rate over angles: 

£=1.57x1 QT 3 B 2 E 2 = aB 2 E 2 , (4) 

where a = b< sin 2 9 >= (2/3)b. Then the electron distribution will remain isotropic, and will simply cut off at £' max (f); a source 
synchrotron spectrum would then cut off exponentially above v(E max (t). 

However, the result we all remember from graduate school is neither of these. A source that turns on at t = with continuous, 
uniformly spatially distributed injection of a power-law distribution of electrons q(E) = JqE~ s electrons erg" 1 s -1 cm" 3 , develops 
a break at energy Et,= 1 / (aB 2 t), where the electron spectrum steepens by one power in s. This corresponds to a steepening of 
the synchrotron spectrum by one-half power in a at vi, = c\B~ 3 t~ 2 with c\ = 1.12 x 10 24 cgs. This relation is frequently used to 
deduce a magnetic-field strength in a synchrotron source of known age. 

The result that synchrotron losses (or inverse-Compton losses, which have the same dependence on electron energy) result in 
the steepening of the particle spectrum by one power and the steepening of the emitted synchrotron spectrum by a half-power, is 
a widely held idea. It is this application that will be generalized below. It is important to remember that the standard derivation 
assumes a distributed injection of electrons in a homogeneous source. For other conditions, it is not correct, as we shall see. 

2. BASIC CALCULATION 

Energy-loss breaks from a synchrotron source in which relativistic electrons are advected systematically away from an injection 
region (such as a wind-termination shock) can be described as being due to shrinking of the effective source size with frequency. 
At a high enough observing frequency, electrons drop below the energy required to emit at that frequency before they reach the 
edge of the object, hence limiting its effective size at that frequency. Thus all results depend on the critical energy an electron 
may have after suffering synchrotron losses and convecting at speed v. Equation H] gives the synchrotron loss rate from a single 
electron. In a constant-density flow (i.e., neglecting adiabatic losses), but allowing the possibility of spatially varying B, we 
generalize the homogeneous results slightly to obtain 
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FIG . 1 . — Schematic of flow geometry. The flow occurs in a tube whose cross-sectional width w grows as a power e of (normalized) distance / from the injection 
radius rg (I = r/rg). 

for the energy an initially infinitely energetic electron would have after moving at v through a field B for a distance r. So any 
injected electron distribution must cut off at this energy. For a one-dimensional flow of electrons streaming at constant speed v in 
a constant magnetic field B, the effective length r(E) of the source is found from 
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An electron of this energy radiates chiefly at frequency v = c m E B 
1970). Then at any frequency v the source has an effective length 
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For synchrotron emission with a spectral index a, the observed flux density would then vary as 

S v oc v~ a r (v) oc u~°~ 1/2 (8) 

- the famous steepening by one-half power in the spectral index. (Note that at no position r in the source is there a break in the 
electron energy distribution N(r,E) to a new power-law, although such a distribution does result after integrating over the entire 
source volume to obtain N(E).) 

This kind of argument can be generalized to inhomogeneous sources, for which the steepening may be larger or smaller than 
0.5. For instance, if the synchrotron emissivity increases with distance from the center, then as the effective source size shrinks, 
more emission will be lost than if the source were homogeneous, and the steepening can be greater than 0.5. 

Here we consider a simple model in which electrons are injected at some initial radius ro in a "jet" whose full width w rises 
with dimensionless distance / = r/ro as Z e : w = wo/ e . A conical jet (or piece of spherical outflow) then has e = 1; a confined jet 
has e < 1 . Then the cross-sectional area increases as Aj_ oc Z 2e (see Figure[T). We shall assume ad-hoc power-law dependences of 
quantities on dimensionless length /, ignoring any transverse variations - so a one-dimensional problem. Let the /-dependence of 
basic quantities be given by 



B = BqI" 
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P = P()1 
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While this is completely genera l, physical const raints may couple the m's. For instance, in the absence of some mechanism such 
as mass loading or entrainment dLvutikovl20 03) to alter the effective density p, conservation of mass gives 

Mass conservation pvA± = const =>■ m p + m v = — 2e. (10) 

In the absence of turbulent amplification or reconnection, magnetic flux will be conserved, giving different relations for compo- 
nents of B parallel and perpendicular to the jet axis: 

Longitudinal (radial) field Bn : BA± = const =>■ mb = -2e (11) 
Transverse (toroidal) field B± : Bwv = const ^>nib= —m v — e = m p + e (12) 

where the last form for transverse B involves invoking mass conservation. Of course, unless the magnetic field is exactly radial, 
any toroidal component will eventually dominate, barring extremely peculiar and probably unphysical behaviors (e.g., m p < -3). 

In this formalism, Equation [5] implies E c oc /-< 1+2m 6-'"i) i n the absence of other energy-loss mec hanisms such as adiabatic 
expansion losse s (non-constant density). In the presence of adiabatic losses, it can be shown (e.g., Ken nel & Coronitil ll984a; 
[Reynolds 1998) that the critical energy is given by 



E c = p 
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(13) 



= -%(1 +2m b -m v + mJ3)l- (l+2 ""--" h> =A E l"' E (14) 

where we have assumed the source is long enough that 1(E) ^> 1. Note that the effects of adiabatic losses have canceled out 
(except for a small change in the integration constant), leaving E c with the same /-dependence as in the constant-density case. We 
have also made the assumption that the integral in Equation Qj] is dominated by the upper limit at /, demanding that m c = 2m b - 
m v +m p /3 >— 1. Otherwise, E c oc / m p/ 3 and radiative losses play no role in the spectral behavior, so that the calculation is not self- 
consistent. Even if this condition is met, we still wish to exclude situations in which the gradients conspire to arrange adiabatic 
gains of electrons as they convect, i.e., > 0. This situation appears both unphysical and unlikely. Note that m c = m p /3- ni£ - 1 . 
These two conditions rule out some volume in the parameter space of (e,m;) and must be checked for any particular choices of 
those parameters. The constraints are related; m c > -1 requires niE < m p /3, so the we ultimately require niE < min(0,m p /3). 
Electrons with energy E c (at position / where the magnetic field strength is B(l)) radiate chiefly at a frequency 

v(E c ) = c m f (1 +2m b -m v + m p /3) 2 n~ 2a+2m *- m * ) = AJ m » (15) 

definingAy andm„. Note m v = m b +2m E = -2-3m b + 2m v . Then 

l(v)=A- 1 }l m »v l l m " =A lV x ' m ". (16) 

We are focusing on conditions such that the source size shrinks with increasing frequency, i.e., m v < 0. It can easily be shown 
that the condition on m c above is equivalent to m v < (m p +m v )- 1. If mass conservation is assumed, m p + m v = —2e and m v < 
always. Otherwise, this condition must be checked, but for reasonable values of the parameters (such as those in all examples 
described below) it is always fulfilled. 

Now to discuss the synchrotron flux, we assume a power-law electron distribution N{E) = KE~ S between energies E\ and 
Eh 3> Ei (we take s > 1). As the flow expands, conservation of electron number gives 

n e = —E\- s txp^Koc pEr 1 oc p 1 ^ 3 = p (s+2) ' 3 = p {2a+3) ' 3 (17) 
s— 1 

since adiabatic losses give individual-particle energies varying as £ oc p x l 3 , at least near Ei, an energy we assume to be too low 
for radiative losses to be important. We have taken s enough larger than 1 that E}~ s 3> Ej^ s . Then we can write the synchrotron 
emissivity (following Pacholczyk 1970) as 

j„ = CjKB l+a V- a = Cj .^ Bl+«/(2a+3)m p /3 + (l + a)« v -a = Ajl m 'V~ a . (18) 

Here Cj(s) = cs(s)(2c[) a in the notation of Pacholczyk; for .5=1 .5, cj = 1 .34 x 10~ 18 cgs, and we have defined another important 
index, mj = (2a + 3)m p /3 + (1 +a)mi,. Assume for the time being that we view the jet directly perpendicular to the axis. Then 
the line-of-sight path length through the jet at any position I is just w(l) = wol e (through the center; averaged over lines of sight 
intersecting a circular cross-section, we obtain the mean line-of-sight path of (n/4)w. We recall that we are assuming all jet 
quantities to be constant in cross-section, that is, along these lines of sight. So if the source is at distance D, the integrated flux 
density S u is given by 

/l f l " f W 7T 7rA ' f l " 

i u dn=- I J^ dAj^ - ju ds=^ i J^ wrQ di( w r>)v- a (19) 

= t rQW 2 dl = * A J r °< 

4D 2 J l 4D 2 {l + 2e + mj ) V ; 

_ TTCjKoB^rowl A (l+2e+(2a+3)m p /3+(\+ a )m b ) v - a -A (2 \) 

4D 2 [l+2e + (2a + 3)m p /3 + (l+a)m b ] ' 
where the last equation defines A, the amount of spectral steepening: 

^ \+2e + nij _ \+2e+m.j _ \+2a + (2a + 3)m p /3 + {\ + a)m b 
m b + 2niE \m u \ 2 + 3m b — 2m v 

We have assumed that m v < 0, and that the flux integral depends on the outer, not the inner, limit of integration, that is, that 

l+2e+m,->0. (23) 

This latter condition can be restrictive. 

While we have assumed a jet seen from the side, since the emission is optically thin the observed flux density should be 
independent of viewing angle. This can be shown explicitly for the case of jets seen exactly end-on, for which the flux integral 
Equation[T9lbecomes 

S » = 4 / ^ ^ wdw I CjK B' +a r oV - a l m >dl. (24) 
u Jo 1 JIM) 

Here the line-of-sight is parallel to the jet axis, beginning at a value of / = dependent on w (if w < H>o,/j = 1), The upper limit 
vfmax is just the value of w at which /, = l v , i.e., w max = w{l(v)) = [l(v)Y. Now we need a slightly more restrictive condition for 



the emission to be dominated by the outer limit of integration l(v): 1 +mj > 0. If this is the case, the two integrals in Equationl24l 
decouple: 

S » = Tk c j K o B o +a rov- a ^ T~ — l+mi ■ (25) 
1U A 1 1 + nij 

But w max = [l(v)Y, so 

S v CX [l(v)] l+2e+mj V~ a CX v ^+2e+mj)/m v ya (26) 

just as in Equationl20l 

These power-law spectra can hold only over a frequency range related to the size range of the source by l(v). For instance, 
for conical, constant-velocity, mass-conserving flow with tangential magnetic field, we have e = 1, m p = -2, and = — 1, giving 
l{v) cx v~ l . Thus a source showing the corresponding value of A (in this case, A = la/3) between frequencies v\ and must 
shrink over a range of sizes given by r\ /rj = vi/v\- In general, sources whose spectra are set by effective variations of size with 
frequency are predicted to have sizes varying as 

l(v) CX V l l m " EE l , 1 /(-2-3m i ,+2m l ,)_ (2?) 

Equationl22lrelates the size exponent 1 /\m u \ to A: 1 /\m v \ = A/(\+2e+nij). (Thus a source with observed A will have a smaller 
rate of shrinkage with frequency for a larger value of l+2e + mj.) The source subtends a solid angle on the sky of 

An = ^ J' (K \w i e )(r di) = ^ [ /(zy)1+e ] ■ (28) 

If the source is only marginally resolved, one may consider an "average" angular size (6) defined by 



(0) ee (Aft)'/ 2 = V^_^[/ M (^>/ 2 (29) 
D Vl + e 

so that {6) cx I /( 1+£ )/ 2m " . For spherical or conical flows, i.e., e= 1, (6) oc v x l m " as before, but for confined flows (e < 1), the average 
angular size decreases more slowly with frequency. If a confined jet is seen end-on, the size variation with frequency is reduced 
even further, as the apparent diameter is now proportional to w mdx cx [l(v)Y instead of l(v): 

9 cx v e l m " = j/A m '> +2m £) = j//(- 2 - 3m 6+ 2m i0 (30) 

This may result in a significantly reduced size effect, since the primary change in emitting volume is a shrinking along the line 
of sight. 

3. SPECIAL CASES 

We can examine a few special cases. First, in the case of a plane constant- velocity flow, we have e = 0,m p = 0, and = 0, and 
we recover A = 1/2. Next, consider tangential magnetic field, and conical outflow (e = 1) with constant density and assuming 
mass and flux conservation. This corresponds to the inner parts of a Kennel & Coroniti (KC) MHD spherical flow. Then m p = 0, 
m v = -2, mi, = m p + e= 1 , and we have 

A — 

which was derived in Kennel & Coroniti (1984b, eq. 4.11b). This already indicates that values of A ^0.5 can be obtained from 
reasonable situations. It also suggests that obtaining values very different from 0.5 may be difficult. The size effect predicted 
above is quite weak: m u = -2-3nib + 2m v = -9, so l(v) cx as shown in Kennel & Coroniti (1984b, eq. 4.10b). This weak 

effect accounts for the very small deviation of A from the homogeneous value of 0.5. 

Now KC models can be divided into two regions: an inner one as above, and an outer one, a constant-velocity flow with p cx r~ 2 
(that is, m v = and m p = -2) and tangential field (nib = m p + 1 = -1). Mass conservation is assumed (m v = -2e-m p ). At lower 
frequencies, the burnoff radius moves in through the outer region, giving 

Aouter = ■ (32) 

However, this is not realized actually, as neither condition m E < or 1 +2e + m ; > is met. First, m E = -(b + 2m\ J + m p ) = +1, 
indicating that E c rises with I. Further, the mj condition gives 3 -2(2a + 3)/3- 1 - a = -la/3, so the flux is dominated by the 
inner parts of the nebula. At high enough frequencies or photon energies that the burnoff radius is at the transition point and 
moves into the v cx r~ 2 region, we obtain the above result 

A ianer = -^. (33) 

Here, 1 +2e+mj ■■ = 4+ a > 0, so the consistency condition is met - the flux is dominated by regions near l(v). For KC's model of 
the Crab, a = 0.6, so Ai nne , =0.51. 

It is straightforward to show directly that spherically symmetric sources obey the same scaling laws with e = 1 - that is, the 
assumptions of thin jets perpendicular to the line of sight still give the correct scaling for spheres. Let the dimensionless radius 
be R ee r/ ro, with injection at R = 1 and burnoff at 

R(V) = A R V l l m » EE A Ri / 1 /(- 2 -3m i ,+2m,.)) (34) 
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- that is, the same expression as for (Sphericity shouldn't change the expression, only, perhaps, the values of the m's.) 
Similarly, we should have the same dependence of j v {R) in the spherical case as we had for j v (l) in the jet case: 

JAR) = CjK {) Bl +a V - a R m ' (35) 

with 

2a + 3 

rrij= — - — m p + (\ + a)m b . (36) 
Then the total flux from this spherically symmetric, optically thin source is just 

4tt f R ^ 47r \R(iA\ m ' +3 

Su = -&\ x ^jAR)R 2 dR = (4ncjK B l +a v- a ) ^.j 3 (37) 

which gives 

A _ m ; + 3 _ i (2a + 3)OT p /3 + (l+a)OT fa +3 
m v 2 + 3mt, — 2m v 

This is the same result as can be obtained by setting e = 1 in the previous expression for A, Equationl22l 

4. BREAKS GREATER THAN 0.5 

Obtaining values of A considerably greater than 0.5, as seems to be required by observations of most PWNe, requires relaxing 
some of the assumptions. The relation from mass conservation of m p = —2e—m v might not hold if some form of mass-loading 
occurs, for instance by evaporation of material from thermal filaments, or entrainment of material from a confining medium. Flux- 
freezing for the magnetic field will not hold in the presence of either turbulent amplification of magnetic field or of reconnection. 
As an example, consider (for algebraic simplicity) the case a = 0, but conical, constant-density flow (so m v = -2), still conserving 
mass. Now in this case, mj = nib, and the = — (3 + 2mj), so 

a 3 + m 3 + m . 6A ~ 3 

nit, — 2(3 + Zrrib) b + inib 3A— 1 

Now we can obtain A = 2/3 with the value nib = _ 1 in this case: the magnetic field drops as the first power of distance (faster 
than if frozen-in and tangential, slower than if longitudinal). This does not seem like an unreasonable possibility. Note that the 
conditions are met: 1 +2e + nij = 3 + m/, = 2 > 0, m E = -1 < 0, and m c = 2m^— m v +m„/3 = > -1. The source size would obey 

oc u 1 '^^ = vVH-W) = V -W . (40) 

To obtain weaker size effects, one is driven to larger values of 1 + 2e + nij, for the same observed A. 

For practical use, it is convenient to consider nib a dependent variable, in terms of the other quantities. First, solve the equation 
for A for nib, in terms of A, m p , e, and a, with no presumed relation between nib and m p . The result is 

A(2m v - 2) + (2a + 3)m „/3 + 1 + 2e 

m b = r-r— : . (41) 

3A- I — a 

Two of the conditions are satisfied if ni£ < min(0,m p /3), or 

1 +2nib-m v > max(0,-m p /3). (42) 

Assuming mass conservation, this becomes 

1 +2e + 2mb+m p > max(0,-m p /3). (43) 

In addition, we still require 1 +2e + m i > 0. 

An application to a particular source (i.e., an object of known a and A, presumably) can then be made by inserting those 
values. Various possibilities for e and m p can be tried; for each value of nib obtained in this way, the condition on /«£ must be 
checked by hand. For exampl e, consider B05 40-693, with a radio spectrum of a = 0.25 (Manchester et al. 1993), and a break 
with A ~ 1 at around 20 fj,m d Williams et al.ll2008l) . A constant-density spherical (or conical) outflow cannot produce this A. 
Inserting the values of a and A into Equation|4Tl and assuming for simplicity e = 1, we find 

l+2m v +1.2m p 

m b = — . (44) 

Then at the price of abandoning mass conservation, we can choose m p = 1 and m v = -2, giving nib = -1-06, or about -1. The 
consistency conditions are all met: m c = 1/3 > — 1, ni£ = — (1 +2m,b— m v ) = -1 < 0„ and l+2e+mj = 2.9 > 0. The source effective 
radius decreases as 

r oc v x l^ lm ^ = i/" 33 (45) 

which may be a serious problem, since we need the slope of -a - A = -1 .2 to hold from about 20 /im to somewhere in the blue 
or near UV - say 0.2 /im, requiring that the source shrink between these two wavelengths by a factor of 4.6 - perhaps unlikely. 
(Though what is shrinking is really the region containing the dominant flux; there could be a faint halo contributing a small 
amount of flux in which a brighter, shrinking core is embedded). A numerical calculation of this model is shown in Figure [2] 
along with observations. The physics which could cause these values of m p and nib is, of course, completely unknown. 
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FIG. 2.— Left: Electron distribution N(EJ) at five positions in the flow model for the PWN B0540-693: / = r/r = 10,30,50,70, and 90. Right: Model 
spectral-energy distribution and observations for B0540-693, reproduced from Williams et al. 2008. Radio: Manchester et al. 1993. IR: Williams et al. 2008. 
Optical: Serafimovich et al. 2004. X-ray: Kaaret et al. 2001. 



m b vs. A: a = 0, e = 1 




m h vs. A: a = 0.3, e = 1 





0.7 

A 



FIG. 3. — Left: magnetic-field index vs. A for a = 0, e = 1, for several values of m p . Mass conservation is assumed. Right: Same, for a = 0.3. Only values of 
A less than the square symbols on each curve satisfy the consistency requirement. 




Figures [3] and |4] plot mi, vs. A for Equation [391 and its generalizations to the pairs (a,e) = (0.3,1), (0,0.5), and (0.3,0.5), 
respectively. Mass conservation is still invoked. The consistency condition on m c places upper limits on Delta (lower limits on 
mt,) shown as the squares on curves on each plot. (The condition 1 + 2e + mj > is met for all curves shown.) It is difficult to 
obtain values of A > 0.7; flows with rapidly dropping density (such as m p = — 2e, for constant-velocity mass-conserving flows) 
seem unable to do so. Substantial deceleration seems to be required, as well as rapid decreases in the magnetic-field strength. 
While there is some parameter space available for accomplishing this, especially for sources with very flat radio spectra, the most 
physically reasonable way to bring about the required deceleration seems to be mass-loading, which also considerably expands 
the available parameter space of source gradients. 
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5. INFERRING PHYSICAL PARAMETERS 

Energy-loss spectral breaks are commonly used to infer source magnetic-field strengths in pulsar-wind nebulae. One requires 
a source age t ; for PWNe in SNR shells, one can use modeling of the shell emission to estimate an age, while in some cases, a 
source size divided by a mean flow speed (estimated one way or another) can give an age estimate. Then one simply assumes a 
homogeneous source for which E c = (aB 2 t)~ l = 637 jB 2 t and from Equation[5] 

B '' = (S) 1/3 ^ 1/3r2/3=0 - 90 ^V 2/3 G. (46) 

(where we have averaged over pitch angles). It is of interest to compare this to the magnetic field that would be inferred for a 
flow model assuming (incorrectly) that the source is homogeneous. Let the source have a size L and break frequency V[ 7 (so that 
L = rQl{vb)). For a flow model, we can deduce the initial magnetic field Bo from EquationfTBlabove: 

L\ m " a 2 1 , 

where /= 1 + 2mb-m v + m p /3. This implies 



Then 



Bo \ V W V r 



Of course, with substantial magnetic-field gradients, B^/Bo can range widely either below or above 1. Some source properties, 
such as the initial ratio of energy input in magnetic field to that in particles (KC's a parameter), require knowledge of Bq. For 
those properties, estimation of source magnetic field from the homogeneous assumption can lead to significant error. However, 
it is possible to show that B/, does give a good approximation to the mean magnetic field averaged over the lifetime of a particle 
moving with the flow, as of course it must. The total magnetic -field energy in a flow model is 

U B =±- [ L/ *r dln(^) 2 Bll** (50) 



871-7, V 2 



1 

32 



w 2 B 2 r 



\+2c+2m b 



1 + 2e + 2nib 

where we have asumed 1 +2e + 2m/, > and L 3> ro. The total volume in the flow is 



(51) 



r.r«(^) 1 »*-M±) M '. (TO 



Then the mean magnetic energy density (ub) is 



2 ) 4(l+2e)Vr 



U B B 2 (l+2e) f L^ 2m " 



<«*> S # =8^172^ U 1 (53) 



and the homogeneous energy density ^(hom) = BjJ&ir satisfies 

«n(hom) _4 /3 (ro\ 4/3 ( l+2e + 2m b \ ( L\ 4(1 -"" )/3 

{ub) J WJ V l+2e JVoJ ( ' 

where the exponent of Lj r ( > has been rewritten using m v = 2m v - 2 - 3m^. 

If the value of t used in the homogeneous relation Equation|46]is the actual transit time of an electron from ro to L, one obtains 
a similar result. That t is given by 

_ fr dl r /• L/ '°,_ m „„ r 1 ( L\ l ~ m < 

''trans — / 

Then 



r"'di = -- - . (55) 

v v J i v 1 -m v Vr 0/ 



4/3 , T x 4(l-m„)/3 
^0 \ [ \ 



Vofeans / V r 



= (l-m l ,) 4 / 3 (56) 



independent of physical parameters. (We are assuming m v < 0, i.e., we exclude accelerating flows.) This means that all factors 
in Equation [54] are of order unity, so that there is not a large discrepancy between the true mean magnetic-field energy density 
and that inferred under the assumption that the source is homogeneous. However, if the source lifetime is used for t (which may 
differ substantially from ? tra nsX or an estimate of f tra ns is made from the measured expansion velocity of the outer boundary of the 
PWN, serious errors may be made in inferring B. 



PWN Synchrotron Spectrum 



PWN Synchrotron Spectrum 
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Log Frequency (Hz) Log Frequency (Hz) 

FIG. 5. — Left: Integrated spectrum for Kes 75 model. Right: Same for MSH 15-52. Flux scales are arbitrary. 

6. NUMERICAL CALCULATIONS AND APPLICATIONS TO OBSERVED SOURCES 

These results can easily be confirmed by numerical integration of the appropriate equations. In particular, Equation [3] can 
be used to find the detailed particle distribution, accounting for the pileup of particles at energies just below E max (t) as they 
migrate down in energy (for s < 2). This pileup can produce a detectable "bump" in the spectrum just below v\,. The numerical 
calculations can also show how sharp a break can be achieved in practice. 

I illustrate these effects with several models. The example parameters for B0540-693 mentioned above (a = 0.25, e = l,m p = 
l,Wi, = — l,w v = -2) produce distribution functions N(E,l) at various points in the flow shown in Figure|2] at positions / = r/ro = 
10,30,50,70, and 90. (The energies are in units of a fiducial energy Ej = aB^o/Vo, the energy an initially infinitely energetic 
electron would have after radiating for a time ro/vo in a magnetic field Bq.) The sharp cutoff energy E c , decreasing down the 
flow, is apparent, as is the spike just below it of electrons formerly above E c . The rising density produces adiabatic gains in 
the density of electrons of too low energy to be subject to radiative losses. Spatial integrati on over these electr on distributions 
produces the model spectral-energy distribution also shown in Figure |2] reproduced from Williams et al. (2008), which fits the 
data surprisingly well, apart from the anomalous X-ray flux. (A technical problem, "pileup" in the Chandra detectors due to 
the bright X-ray pulsar in B0540-693, makes the absolute determination of the X-ray flux of the nebula difficult; see Petre 
et al. 2007.) For the relatively steep low-frequency spectrum of B0540-693 (a = 0.25), the "bump" from integrating over the 
spikes of Figure [2] i s barely noticeable, but it is much more obvious for a flatter input spectrum. The model assumes a source 
radius L ~ 1.3 pcdWilliams et alJl 2008), but ro is a free parameter. If rn is the pulsar wind termination shock, we might expect 
vo = c/3 (Kennel & Coroniti 1984a). The break frequency (Equation [15]) constrains the remaining combination r^B^. The model 
of Figure [2] takes ro = L/100 and Bo = 2.4 x 10~ 3 G. The radio flux (Equation |2TI) then sets the combination WqKo; the model 
takes wo = ro/ 10 and Kq = 2. 1 x 10~ 8 cm" 3 erg 5 . The break frequency calculated from Equation[l5]is about 6 x 10 12 Hz, about a 
factor of 5 lower than the intersection of the extrapolations from low and high frequencies. This is due to the approximation that 
electrons radiate entirely at their peak frequency v m , an approximation not made in the analytic calculations; in general, break 
frequency predictions will be low by a factor of several, depending somewhat on the value of s. 

Chevalier (2005) summarizes spectral indices for several PWNe, including Kes 75 (a = 0=>s = 0, A = 1) and MSH 15-52 
(a = 0.2 =>■ s = 1.4, A = 0.85). Such large values of A typically require relaxing either mass or flux conservation, or both. For 
Kes 75, the values e = 1, m p = 1, m v = -2, and nth = _ 1 (the same as for the B0540-693 model, except a = 0) predict A = 1.0. 
The consistency condition \ +2mi,-m v > max(0,-m p /3) is met. Figure [5] (left) illustrates the integrated spectrum. The "bump" 
is quite prominent; the flux at the peak around 10 1 Hz is 4.4 times that at 1 GHz. The slope above the break is 0.96 between 
0.4 and 4 keV, close to the analytic value of 1.0. Equation[27]gives the frequency-dependence of the source size as Z max oc i/" 1 / 3 , 
sufficiently slow that it might be hard to detect. For MSH 15-52, Figure (right) shows a calculation for s= 1.4, e= 1, nt p = 1, 
m v = -2.22, and mi, = — 1, predicting A = 0.85. The consistency condition is again met. The "bump" is still perceptible. The 
predicted value of A is reproduced exactly. The size effect is even slighter: Z max oc v 29 . A factor of 1 1 frequency range would 
be required to see the source shrink by a factor of 2. For an actual well-resolved source, the "size" would need t o be measured with 
some r elatively coarse quantity such as the 50% enclosed power radius or FWHM, as cited for the Crab by iKennel & Coronitil 
(Il984bl) . 

7. SUMMARY OF RESULTS 

Here I collect the principal results and consistency requirements. The basic result is the expression for a^-ai = A, the amount 
of spectral steepening, Equation [22l 

^ _ 1 + 2e + nij _ 1 + 2e + (2a + 3)m p /3 + ( 1 + a)m b 
nii, + 2m E 2 + 3wfc — 2m,, 

This expression holds if several consistency requirements are met: Me < min(0,m p /3), where me = -1 - 2mt, + m v , so that the 
burnoff energy at position / depends on /, and E c drops with I; and 1 + 2e+ (2a + 3)m p /3 + (1 + a)nib > 0, so that the integrated 
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flux density S„ depends on the outer limit of integration. Finally, the effective source size should shrink with frequency: m v < 0, 
a condition always met in the presence of mass conservation, and almost always met for reasonable parameters otherwise. If the 
conditions are met, the source size (some measure of the region from which the bulk of the emission originates) decreases with 
frequency as 

l(i/)ai/ 1/ " , "Ei/- 1 / (2+3 ""- 2 '"' ) . (58) 
if seen more or less from the side; if the flow is nearly along the line of sight, the expression becomes 

6 oc i/ /m " = ^A***-^) (59) 
(which is the same for spherical or conical flows where e = 1). 

8. CONCLUSIONS 

My basic conclusion is just that synchrotron-loss spectral breaks differing from 0.5 can be produced naturally in inhomo- 
geneous sources. I have treated the inhomogeneities resulting from flows, which seem most natural, using simple power-law 
parameterizations, but more complex functional dependencies can be treated the same way. Other types of inhomogeneities 
may be possible as well. These results are most straightforwardly applied to PWNe or knots in extragalactic jets, but may have 
applications wherever bulk flows of relativistic material are involved. In particular, energy-loss breaks seen in gamma-ray burst 
afterglows (e.g., Sari et al. 1998; Galama et al. 1998; and much later work) may provide opportunities for the application of 
these results. For nearby sources, the simplest test of the models is the detection of the size effect; every model predicts both a 
particular A and some rate of decrease of size with frequency (really of volume, since the size decrease may take place along the 
line of sight). For the same A, some variation in the strength of the size effect is possible. 

Fortunately, assuming that an inhomogeneous source is actually homogeneous does not drastically alter the inferred mean 
magnetic field (averaged over the history of a fluid element), but since there are by assumption large gradients of most quantities, 
local values of the magnetic-field strength, such as those at the injection radius, may depart substantially from the mean values. 
This may affect inferences of the KC magnetization parameter a. Furthermore, the inference requires knowledge of the actual 
flow time across the source - knowledge that may be hard to come by in the presence of large velocity gradients. 

I gratefully acknowledge the hospitality of the Arcetri Observatory of the University of Florence, where this work was begun. 
This work was also supported by NASA through Spitzer Guest Observer grants RSA 1264893 and RSA 1276758. 
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